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The Anderson transition in three dimensions in a randomly 
varying magnetic flux is investigated in detail by means of the 
transfer matrix method with high accuracy. Both, systems 
with and without an additional random scalar potential are 
considered. We find a critical exponent of v — 1.45 ± 0.09 
with random scalar potential. Without it, v is smaller but 
increases with the system size and extrapolates within the 
error bars to a value close to the above. The present results 
support the conventional classification of universality classes 
due to symmetry. 
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Since the pioneering work of Anderson jjj, the 
disorder-induced metal-insulator transition, which is one 
of the most fundamental quantum phase transitions in 
condensed matter physics, has attracted considerable at- 
tention Depending on the symmetry, the critical 
behavior of this Anderson transition(AT) is convention- 
ally classified into three universality classes: the orthog- 
onal, the unitary and the symplectic class Q. Systems 
invariant under spin rotation as well as under time rever- 
sal belong to the orthogonal class. Unitary systems are 
characterized by the absence of time-reversal symmetry, 
due to, for instance, a magnetic field. Systems without 
spin rotation invariance belong to the symplectic class. 

The AT in a homogeneous magnetic field has been 
studied extensively for many years, mainly in connection 
with the quantum Hall effect ||. In two dimensions(2D) 
in the presence of a strong magnetic field, the AT is 
marginal. States at the centers of the Landau bands are 
critical and all the other are localized. At the band cen- 
ters, the localization length diverges with the exponent 
v ~ 2.4 IJ . In 3D, there exist extended states and the 
AT takes place The latter has recently been re- 

analyzed and the critical exponent for the localization 
length has been determined to be 1.43 ± 0.06 flC(| . 

When the magnetic field is uniform in space, random- 
ness is introduced by a random scalar potential. On the 
other hand, in recent years, there has also been consider- 



able interest in 2D systems subject to a spatially random 
magnetic field, mainly in connection with the fractional 
quantum Hall effect |11|. The random magnetic field in- 
troduces randomness as well as the absence of invariance 
under time reversal in a system. 

In 3D, the AT in the presence of a random vector po- 
tential and without a random scalar potential, has been 
investigated numerically. The data suggested that 
the mobility edge is very close to the band edge. The 
exponent for the localization length has been estimated 
to be v ss 1 1 12 1 which is considerably smaller than that 
in the case with an additional random scalar potential 
and in a uniform magnetic field. It has also been re- 
ported |||n| that in the presence of a random scalar po- 
tential, the critical exponent has a universal value, irre- 
spective of whether the magnetic field is uniform or ran- 
dom. This seemed to indicate that the AT in a random 
vector potential but without random scalar potential is 
different from the one with a random scalar potential. 
It should be noted that this would question the valid- 
ity of the above conventional classification of the AT be- 
cause in both cases the time reversal symmetry is broken 
and hence these two systems should belong to the same, 
namely the unitary universality class. 

The critical exponent v w 1 for the 3D system with a 
random magnetic field has been obtained by the finite- 
size scaling method jl4). This method has been applied 
successfully to analyze the critical behavior of the AT || . 
In most cases, however, the numerical analyses have been 
restricted to energies near the band center. It has been 
reported that systematic scaling behavior has not been 
clearly observed for energies away from the band center 
|l5| , ^6| . For the 3D system with a random magnetic field, 
the mobility edge lies quite close to the effective band 
edge [^2|. It is therefore imperative to investigate the 
present problem with considerably higher accuracy and 
to examine carefully whether or not the scaling behavior 
is modified by adding a random scalar potential. 

Recently, high-accuracy scaling analyses of the Ander- 
son transition have been performed by several authors 



1 



0,0. It has been concluded |L0l that, by reducing the 
errors of raw data to 0.1 ~ 0.2%, one can numerically 
distinguish the unitary from the orthogonal class, which 
was impossible when only low accuracy data with *~ 1% 
accuracy were used. 

Encouraged by this recent success, we have started a 
numerical high-precision finite size scaling project, in or- 
der to clarify the above mentioned discrepancy between 
the critical exponent of the AT far away from the band 
center induced solely by randomness in a vector potential 
and the exponent obtained for other unitary systems. 

We found a clear systematic dependence of the expo- 
nent on the system-size in the former case which would 
introduce corrections to scaling. These would become 
smaller than the statistical error only for system sizes 
larger than those which are presently achievable. We esti- 
mate the asymptotic limit for the exponent to be v ~ 1.4, 
when the transition is near the band edge. This behavior 
is not changed significantly when shifting the mobility 
edge by adding a weak random scalar potential. 

The model is defined by the Hamiltonian Il2) 

<i,j> i 

where Cj(Ci) denotes the creation(annihilation) operator 
of an electron at the site i of a 3D cubic lattice. Ener- 
gies {Vi} denote the random scalar potential distributed 
independently and uniformly in the range [-W/2, W/2]. 
The Peierls phase factors exp(i#jj-) describe a random 
vector potential or magnetic field. We confine ourselves 
to phases {0ij} which are distributed independently and 
uniformly in [— 7r, 7r]. The hopping amplitude t is as- 
sumed to be the energy unit, t = 1. 

We consider quasi- ID systems with cross section M x 
M Tnc Schrodinger equation Hip = Etp in such 

a bar-shaped system can be rewritten using transfer ma- 
trices T„(2M 2 x 2M 2 ) 

/ Vn+l \ _rp f Ipn \ rp _ ( E-H n -I\ 

\ ^ J-^y^n-xJ' '"-{I ) 

(2) 

(n = 1,2,...) where ip n and H n denote the set of coef- 
ficients of the state t/j and the Hamiltonian of the n— th 
slice, respectively. The identity matrix is denoted by I. 
The off-diagonal parts of the transfer matrix T n can be 
expressed by the identity matrix because the phases in 
the transfer-direction can be removed by a gauge trans- 
formation |L2| . The logarithms of the eigenvalues of the 
limiting matrix T 

n n 

T = n 1 ™[(Il^) t (Il^)] 1/2 " ( 3 ) 

i=l »=1 

are called the Lyapunov exponents. The smallest Lya- 
punov exponent Am along the bar is estimated by a 



technique which uses the product of these transfer matri- 
ces §Q The relative accuracies for the smallest Lya- 
punov exponents achieved here are 0.2% for M < 10 and 
0.25% ~ 0.3% for M = 12. The localization length £ M 
along the bar is given by the inverse of the smallest Lya- 
punov exponent, £m = 1/Am- 

The assumption of one-parameter scaling for the renor- 
malized localization length Am = £m/M implies 

Am = /(£/M), (4) 

where £ = $(E, W) is the relevant length scale in the limit 
M -> oo flj. Near the mobility edge E C (W), £ diverges 
with an exponent v as £ ~ x~ u with x = (E — E c ) / E c . If 
the transition is driven by the disorder W at a constant 
energy, x = (W c - W)/W c . At the mobility edge, Am 
becomes scale-invariant. The quantity Am is a smooth 
function of E and W , and we can expand it as a function 
of x as 

oo 

A M =A C + J2MM 1/U x) n . (5) 
n=l 

By fitting our data to the above function, we can de- 
termine the critical exponent v and the mobility edge 
accurately. In practice, we truncated the series (^|) at 
the third order. 

We used the standard x 2 -fitting procedure |Q . In or- 
der to check the goodness of the fit, we also evaluated the 
probability Q that the \ 2 wm exceed the minimum value 
Xmin actually obtained by the fit. The probability Q is 
evaluated via the incomplete gamma functions and the 
condition Q > 0.001 is often regarded as an acceptable 
condition for the fitting function If the value of Q 
is too small, in other words, the minimum value of the 
X 2 is considerably large, there may be systematic devia- 
tions of the numerical data from the fitting function. In 
the recent work fic[l on the Anderson transition in 3D or- 
thogonal and unitary systems, it has been demonstrated 
that the above fitting function up to the third order is 
in fact valid. The error bars are estimated by using the 
Hessian matrix and the confidence interval is chosen to 
be 95.4%. 

We consider first the AT at the band center E = 
in the presence of a strong random scalar potential as 
well as random vector potential. The renormalized lo- 
calization lengths Am evaluated for the disorder W in 
the range [17.8, 19.8] and sizes M = 6, 8, 10 and 12 are 
shown in figure 1. The above described fit yields the crit- 
ical exponent v = 1.45 ± 0.09 and the critical disorder 
W c = 18.80 ±0.04. The renormalized localization length 
A c at the critical point is 0.558 ± 0.003. The value of Q 
is ~ 0.89 which confirms the validity of the fitting func- 
tion (^|) and thus of the one-parameter scaling behavior 
of Am in this range of the disorder W. The error bars of 
the present results are at least a factor of 3 smaller than 
those of the previous estimates || . 
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FIG. 1. The renormalized localization length Am as a func- 
tion of disorder W for different sizes. The crosses, triangles, 
the squares and the diamonds correspond to M = 6, 8, 10 
and 12, respectively. Inset: The scaling function. 

Next, we concentrate on the AT in the random mag- 
netic field but without random scalar potential (W = 0). 
In the previous work fl2|| , the energy range used for the 
scaling analysis was assumed to be 4.3 < \E\ < 4.5. This 
is very close to the band edge and thus the density of 
states(DOS) is rapidly decreasing |fL2[ . To get rid of the 
influence by this rapid change of the DOS, the energy 
window for the scaling analysis should be taken to be as 
small as possible. We therefore choose calculated data 
for 4.39 < E < 4.44. The energy window is 4 times 
smaller than that in |12|. The numerical data for Am are 



shown in figure 2. The transition can be located around 
E w 4.415. We then fit the data for M = 6, 8, 10 and 
12 to the function (|). This yields E c 4.414, v « 1.2 
and A c 0.51, which is consistent with the previous es- 
timates |l2] ]. It is found, however, that in the present 
case, the value of Q turns out to be very small , namely 
~ 10~ 14 , in contrast to the above considered case of the 
band center. This shows that systematic deviations of 
the data from the fitting function are very likely to exist. 
We therefore have to analyze the numerical data much 
more carefully. 

In order to get insight into the origin of the deviations 
of the data from the fitting function, we carried out the 
fits using different combinations of system-sizes Mi and 
M 2 = Mi + 2 for Mi = 6, 8 and 10 (table I). Although 
the crossing point is almost size independent, the expo- 
nent shows a systematic dependence on the system cross 
section. For M = 6 and 8, in particular, the exponent 
is close to 1, while for the other two cases it is around 
1.3 which is close to the value estimated for a system in 
a uniform magnetic field. In addition, for M = 6 and 
8 the value of Q becomes much smaller than those for 
the other two cases, which could imply that the fitting 
function is not working well for these smaller sizes. 
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FIG. 2. The renormalized localization length Am as a 
function of energy E for W = 0. The crosses, triangles, 
the squares and the diamonds correspond to M = 6, 8, 10 
and 12, respectively. Inset: The logarithm of the derivative 
Ci ee dA M /dE for W = at E = 4.414 w £ c (nlled triangles) 
and for W = 1 at E = 4.451 « i5 c (filled diamonds) as a func- 
tion of InM. The solid line and the dashed lines represent 
the slope with v — 1.45 and v = 1, respectively. 

We also analyzed the data in a different way. We per- 
formed a third order polynomial fit for each size and es- 
timated the derivative C\ = dh.M / <1E\e=e c at the crit- 
ical energy E c — 4.414. The value of Q in this case is 
larger than 0.1 for any size, which means that the 3rd 
order polynomial fit itself is working fairly well. We plot 
In | Ci | as a function of InM (figure 2, inset). If scaling 
works, the slope is related to the exponent by \ jv. It is 
clear from figure 2 (inset) that the slopes deviate from 1 
and are likely to approach the value v ~ 1.45, when M 
is increased. These two analyses suggest that the present 
system-sizes may not be large enough to observe clear 
one-parameter scaling behavior for W = 0. 

To see whether these features are specific to W = 0, 
we performed the same analysis for a system with weak 
scalar randomness W = 1. The results of the fits using 
different combinations of system-sizes and by the third 
order polynomial fits are listed in table II and shown in 
figure 2 inset, respectively. In table II, we see that the 
exponent for M = 6 and 8 again turns out to be close to 1 
and deviates from those for the other two cases, although 
the crossing points are rather stable with respect to the 
change of sizes. It should also be noted that the value 
of A c is almost the same as for W — 0. In the inset of 
figure 2, we can observe again the deviation of the slopes 
from 1 as the size is increased. Thus, the system with a 
weak random scalar potential shows similar behavior as 
that without a random scalar potential. 

The results for W = 1 and W = indicate that the cor- 
rection to scaling is not specific to the case W = 0. It is 
natural to expect that this is due to the fact that the tran- 
sition is near the band edge. In both cases, W — 1 and 
W = 0, the fits using the size M = 6 give smaller values 
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of the critical exponent v. We conclude that this is the 
reason why the exponent found in the previous work ]T^ ] 
was smaller. The present analysis for W — 1 and W — 
also shows that finite-size corrections to scaling exist for 
the presently achievable system sizes (6 < M < 12), 
which might be the reason of the discrepancy between 
A c near the band edge and that at the band center. 

In summary, we have re-investigated in detail the An- 
derson transition in the random magnetic field with and 
without random scalar potential. We have evaluated the 
localization length along quasi-lD systems with high ac- 
curacy and examined the scaling behavior of the renor- 
malized localization length near the transition. We have 
confirmed the one-parameter scaling behavior for the 
transition at the band center with a relatively strong ran- 
dom scalar potential and found the exponent 1.45 ±0.09. 
This value agrees well with the recent precise results 
for systems in a uniform magnetic field jnj. We have 
also performed the finite-size scaling analysis for both, 
the system without random scalar potential (W = 0) 
and with a weak random scalar potential (W=l). In 
both cases, deviations of numerical data from the scaling 
ansatz are found, especially for smaller sizes. As the size 
is increased, the exponent is more likely to approach to a 
value around 1.4 rather than the value of 1. In particular, 
no evidence for the exponent 1 has been found. 

On the basis of the present results, we conclude that 
there exists no evidence that the critical behavior in a 
3D system in a random magnetic field is different from 
that for other unitary systems. This supports the con- 
ventional classification of the AT by universality classes 
due to symmetry. 
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(M 1 ,M 2 ) Ac v 


E c 


Q 


(6,8) 0.514 ±0.005 1.05 ±0.07 
(8,10) 0.516 ±0.007 1.26 ±0.09 
(10,12) 0.51 ±0.01 1.32 ±0.12 


4.414 ±0.001 
4.414 ± 0.001 
4.414 ± 0.001 


~ 10~ b 
~ 0.89 
~ 0.99 


TABLE I. Results of the fits for different sizes in 
sence of random scalar potential (W = 0). 


the ab- 


(Mi,M 2 ) A c v 


E c 


Q 


(6,8) 0.510 ±0.005 1.09 ±0.08 
(8,10) 0.519 ±0.008 1.36 ±0.12 
(10,12) 0.51 ±0.01 1.34 ±0.14 


4.451 ±0.001 
4.450 ± 0.001 

4.452 ± 0.002 


~ 0.91 
~ 0.56 
~ 0.95 



TABLE II. Results of the fits for different sizes in the pres- 
ence of the weak random scalar potential (W = 1). The 
energy window is taken to be 4.425 < E < 4.475 and the 
same number of data points as the case of W = are used 
for the scaling analysis. 
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